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In previous works, a squared mass operator and a more common centre of 
frame Hamiltonian -ffcM for the quark-antiquark system were obtained by using 
the only assumption that In W {W being the Wilson loop correlator) can be written 
in QCD as the sum of its perturbative expression and an area term. We evaluated 
the spectrum of these operators by neglecting all the spin-dependent terms but 
for a spin-spin term. We succeeded in reproducing rather well the entire meson 
spectrum and for the light-light systems we found straight Regge trajectories with 
the right slope and intercepts. A first attempt to study the contribution of the 
spin dependent terms has been made. 



1 Deduction of the quadratic and linear mass operators 

It was shown in previous works that starting from proper two-particle and 
one-particle gauge invariant Green functions H^^{xi, X2] j/i, 2/2) and H^^{x — y) 
it is possible to obtain a first principle Bethe-Salpeter and Dyson-Schwinger 
equation respectively □. These gauge invariant Green functions are written in 
terms of a "second order" quark propagator A'^(a;,?/;A) in a external gauge 
field A^, defined as follow by the iterated Dirac operator 

{D^D'^ + m'~]^g a^^^F^^Wix, A) = - 5\x - y) (1) 

wither^'' = (z/2)[7^, 7'^]. By a Feynmann-Schwinger representation, A°'(x,?/;A) 
can be written in terms of the Wilson correlator 

W[T]^^-{TvVeMwj>^dx^A^})- (2) 

In order to allow the evaluation of W we have to introduce some assumptions. 
In the modified area law (MAL) model it was assumed 

Zlniy = z(lniy)pert+CT5min (3) 

where the first quantity is supposed to give correctly the short range limit 
and the second the long range one. If one replaces the minimal surface 5*111111 
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by its "equal time straight line approximation" and omits in the path inte- 
gral representation of H^^{xi,X2,yi,y2) and H^^{x — y) the contributions to 
i In coming from gluon lines and fermionic lines involving points of the 
Schwinger strings, by an appropriate recurrence method, an inhomogeneous 
Bethe-Salpeter equation and a Dyson-Schwinger equation can be derived re- 
spectively. The corresponding homogeneous BS-equation can be written as 

, , f d'^u - , 1 k + u 1 k + Us 

Hi{^P + k)a''^p(u)a'H2{-^P + k) (4) 

and in terms of the irreducible self-energy H{k) = Ho{k) + iHo{k)r (k) H (k) 
the DS-equation can be written as 

m = j^Ja,{k-l-^i'-±i)a^H{l)-'. (5) 

The kernel lab is completely known and because of eq. (|^) it is given as the sum 
of a perturbative term and a confinement term. In order to solve eq.(^) we 
should solve eq.(||) and substitute H into eq.(U), then we should solve eq.(0). 
We can simplify the problem by approximating H with a free propagator 
—ijij? — rr?) and performing an instantaneous approximation that consists 

in setting ko = k'^ = (l/2)[r72(u'i + w[) — 771 (u'2 + w'2)] with wj = -I- , 

Wj = 'iJiTij + k'2 , r]j — mj /{nil + "^2) and j — 1, 2. By performing explicitly 
the integration over fcp and fco we obtain 

{Wi+W2f VmsO^) + (6) 
[ d?'k' Iwi +W2 finst/i 1 /N h"l+'^2 a n i\ b 2 n\ 

J (27r)-^ V 2wiW2 y 2wiW2 

with 'y3p(k) ~ yj2wiW2/{wi + W2) J^^dko^p{k). Eq.(||) is the eigenvalue 
equation for the squared mass operator 

m2 ^ + U (7) 
with Mo = i/mfn? + -f and 

(27r)-^ V 2w\W2 y 2wiW2 
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The quadratic form of eq. obviously derives from the second order character 
of the formahsm we have used. In more usual terms one can also write 



HcM = M ^Mo + V- 



(9) 



with 

(kiyik') = 



1 f'^^'rW \c'\ 

-(k|[/|kO = ^ ' / (10) 



(27r)3 A^/ wiW2w[w2 



In eq.(^ the dots stand for higher order terms in as and cr; kinematic factors 
equal to 1 on the energy shell have been neglected. 
From eqs.(H) and one obtains explicitly 



(k|C/|k') = 



' {Wi + W2){'W[ + W2) J 4 Q!s 
WiW2w'iW'2 \ 3 TT^ 
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■[q2o{ai ■ Q) - 910(0:2 • Q)] 



--a,-a2 + -\^-rT,-CT2 

-/(^e^^-^'/""(r,q,?io,92o)} 



2Q2' 

(Q-cri)(Q.cr2)\ , 1 



4Q2 



(ai-Q)(a2-Q) 



(11) 



with 



J'"''*(r,q, 910,920) = 



910 + 920 



920 V 9?o - Qt + 9?o V 920 - q| 



9?o9lo/ • Iqxl , . Iqxl 

--; f-{ arcsm h arcsm 



iqxl 



910 



920 



920 



Vi'io - qx 



(r X q • o-i 



+z9io(r • ai)) + 



910 



:(r X q - (T2 - i92o(r ■ "2)) 



(12) 



V 120 ^ qx 

here a'^ = j^j^, = ^ £*="™[77, 77'], Q = k - k', q = (k + k')/2 and 
9jo = {wj + Wj)/2. 

Due to eq.(p^) the potential V can be obtained from U as given by eqs.(|l 
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( p^ ) simply by the kinematic replacement 
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Figure 1: Heavy- heavy and heavy-light quarkonium spectrum. Circlets are our results by 
the HcM operator, crosses are our results by the operator. 



2 Quarkonium spectrum and Regge trajectories 

We computed the heavy-heavy, heavy-Hght and hght-hght quarkonium spec- 
trum and the the Regge trajectories for the hght-hght quarkonium systems 
by using both the hncar mass operator iJcM and the quadratic mass opera- 
tor Af^. We solved their eigenvalue equation by the Rayleigh-Ritz variational 
method using the three-dimensional harmonic oscillator basis 13. In the linear 
potential V and in the quadratic interaction U we have taken into account 
only the velocity dependent terms in order to compute the baricentral meson 
masses and the spin-spin ai ■ <Ji term in order to compute the hyperfine split- 
ting. We have then neglected all the other spin dependent terms that will be 
studied later. As shown in figsjl]-|| our theoretical meson masses are in good 
agreement with all the experimental data but for the light pseudo-scalar states 
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Figure 2: Light-light quaxkomium spectrum and Regge trajectories. Circlets and dotted 
lines are our results by the i?cM operator, crosses and dashed lines are our results by the 
operator. 
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which are related to the chiral symmetry problematic. We found also straight 
Regge trajectories with the right slope and intercepts. 

The introduction of a running coupling constant was especially required in 
the computation of the eigenvalues. However it was introduced also in 
the HcM case. We chose the simplest perturbative running coupling constant 
as(Q) = (47r)/[(ll - |Nf) In (Q^/A^)], with Nf = 4, A = 200 MeV and a 
horizontal cut at a value as(0) in order to avoid the singular point = A^. 
We chose: the heavy quark masses nih and mc, in order to reproduce exactly 
the masses of J/vf and T; the running coupling constant cut as(0), in order to 
reproduce the charmonium hyperfine splitting; the string tension a, in order 
to obtain a right slope for the uu Regge trajectories. No fit was made on the 
light-heavy quarkonium systems. „ 
The light quark masses were fixed on typical current values mu — — 
10 MeV and mg = 200 MeV both for the first order formalism and for the sec- 
ond order formalism. 

In the first order formalism we adopted the following parameters crEl: mc — 
1.40 GeV, mb = 4.81 GeV, as(0) = 0.363, a = 0.175 GeV^iJvhereas in the sec- 
ond order formalism we adopted the following parameters mc = 1.394 GeV, 
mb = 4.763 GeV, a^iO) = 0.35, a = 0.2 GeV^. 

A first attempt to study the contribution of the spin dependent terms has 
been made. For the heavy-light systems a lot of spin dependent terms are 
negligible and we are left with a small number of them. The problem becomes 
very similar to the study of the hydrogen-like atom spectrum by means of the 
quadratic form of the Dirac equation eI. This part of the work is in progress. 
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